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INTRODUCTION 


In  connection  with  a  r.ew  method  (Reference  1)  for  estimating  the  struc¬ 
tural  loads  applied  to  an  airplane  by  flight  in  turbulence,  the  FAA  has  supported 
this  brief  study  of  the  estimation  of  airplane  stability  derivatives.  Numerical 
values  of  the  latter  are  needed  in  formulae  for  the  loads,  and  the  accuracy  of 
the  loads  estimate  directly  depends  on  the  accuracy  of  the  stability  derivatives. 

The  new  method  for  gust  loads  was  a  simplified  one,  intended  for  use  by 
manufacturers  of  light  general  aviation  aircraft,  and  so  the  stability  derivatives 
of  interest  would  be  those  of  light  planes  with  straight,  high  aspect -ratio  wings 
and  conventional  tails,  and  propeller  propulsion.  The  Navion  airplane,  which 
has  been  extensively  operated  and  tested  at  Princeton  University,  was  to  be 
the  subject  of  the  stability  derivative  study. 

The  derivatives  of  the  Navion  were  tc  be  predicted  by  several  different 
standard,  or  conventional,  methods.  They  would  be  compared  with  full-scale 
wind-tunnel  data  recently  available,  and  finally  with  values  deduced  from  dy¬ 
namic  responses  recorded  in  a  special  series  of  flight  tests.  The  latter,  of 
course,  would  be  the  norm  against  which  the  other  estimates  would  be  evaluated 
for  accuracy. 

The  standard  methods  to  be  considered  were  selected  more -or -less  arbi¬ 
trarily  to  be  commonly  available  textbooks,  NACA  reports,  and  USAF  DATCOM. 
Although  all  of  those  references  call  attention  to  the  necessity  for  judgment  and 
experience  in  applying  their  methods,  formulas  and  charts;  in  the  context  of  the 
present  study  it  was  deemed  best  to  follow  literally,  without  modification,  the 
procedures  presented.  The  estimates  derived  are  what  would  be  calculated  by 
an  inexperienced  engineer  -  one  who  could  follow  instructions  without  mistakes, 
but  one  without  the  background  or  confidence  to  seek  other  data  or  to  make  arbi¬ 
trary  variations,  as  might  seem  desirable  to  an  "old  band"  at  vehicle  aero¬ 
dynamics.  The  resulting  accuracy  is  what  might  be  expected  from  a  new 
engineering  graduate  given  the  task  to  estimate  the  derivative?  without  guid¬ 
ance  from  a  senior  supervisor. 
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The  numerical  values  of  parameters  predicted  by  the  three  textbooks 
were  calculated  first  by  a  group  of  students  and  checked  by  members  of  the 
research  staff  at  Princeton  University.  Methods  and  instructions  were  fol¬ 
lowed  literally.  Where  alternate  methods  were  given,  the  recommended  ones 
were  selected. 

The  USAF  Datcom  estimates  were  prepared  by  Douglas  personnel,  by 
request  of  the  Air  Force.  The  Datcom  methods,  of  course,  cover  a  wide  range 
oi  airplane  configurations  with  little  emphasis,  actually,  on  conventional  light  - 
plane  types.  Some  of  the  methods  are  not  directly  applicable,  so  that  good  ac¬ 
curacy  in  this  reference,  for  this  problem,  is  not  necessarily  to  be  expected. 

The  lateral -directional  parameters  listed  for  NACA  TR  1098  (Refer¬ 
ence  2)  do  not  include  the  important  Cy£,  Cn^ ,  and  Cnr  derivatives.  The 
vertical  tail  contributes  greatly  to  these,  and  although  the  reference  describes 
methods  for  correcting  empirical  data,  the  basic  data  are  not  given  in  the  refer¬ 
ence.  Consistent  with  the  ground-rules  of  this  investigation,  the  necessary 
reference  to  other  material  is  disallowed,  and  no  estimates  are  possible. 

The  wind  -tunnel  estimates  of  parameters  are  data  taken  in  full-scale 
tests  at  Langley  Research  Center  in  1969  (Reference  3).  The  airframe  was 
originally  an  actual  Navion,  modified  for  variable  tail  incidence,  special  flap 
deflections,  and  reverse  thrust  propeller.  For  the  wind  -tunnel  tests,  an 
electric  motor  was  installed.  The  two  conditions  involved  correspond  to  a 
moderate  cruise  at  150  mpL,  and  a  landing  pattern  condition  with  1/2  flap 
and  power  for  level  flight  at  90  mph. 

The  flight  values  of  parameters  are  derived  from  special  flight  tests 
of  the  Princeton  variable -stability  Navion,  N91566.  Techniques  and  special 
considerations  are  discussed  in  subsequent  sections. 

In  the  gust  response  calculations  (Reference  1)  both  longitudinal  and 
lateral -directional  stability  derivatives  are  involved.  The  following  of  the 
two  groups  are  required: 
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Although  they  are  not  involved  in  the  simplified  calculations  of  Reference  1, 
two  other  derivatives  are  of  general  interest  in  stability  calculations,  and  they  are 
included  in  the  following  study,  for  completeness.  They  are 


The  Navion  Airplane 

The  physical  characteristics  of  the  Navion  are  presented  by  the  three - 
view  drawing  of  Figure  1  and  the  following  tabulation  of  dimensions  and  inertial 
properties.  Two  flight  conditions  are  considered,  identified  below  as  I  and  II. 
The  first  is  a  moderate  level-flight  cruise  condition  with  flaps  and  gear  up,  and 
about  75%  power;  whereas  the  second  is  a  landing  pattern  condition  with  half -flap 
deflection  and  power  for  level  flight.  The  general  characteristics  are  as  follows. 
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TABLE  I  -  NAVION  AIRPLANE  DIMENSIONS 


Wing 

Area,  S 

1 84  ft2 

Sweep 

2°  59'46m 

Aspect  Ratio,  A 

6.  04 

Taper  Ratio,  X 

.54 

Mean' Aerodynamic  Chord,  c" 

5.  7  ft 

Dihedral 

7.5° 

Incidence  Root,  iw 
wr 

+2° 

Incidence  tip,  iw 

-1° 

Airfoil  tip 

NACA  6410  R 

root 

NACA  4415  R 

Horizontal  Tail 

Area 

43  ft2 

Sweep 

6° 

Aspect  Ratio 

4.  0 

Taper  Ratio 

.67 

Airfoil 

NACA  0012 

,0 

Incidence 

-3 

Vertical  Tail 

_  5 

12.5  ft2 

NACA  0013.2  Mod  (root)  NACA  0012.04  Mod 
2o  (tip) 


Area  (above  horizontal  stabilizer) 

Airfoil 

Fin  Offset 


5. 


Propeller  Characteristics 

Diameter  84" 

Number  of  blades  2 
Side  force  factor  100 


Power  Plant 

Continental  engine;  Model  #1052 OB 
HP  Rating  285  at  take-off  at  2700  RPM 


Control  Surfaces 

Surface 

Area  (ft2  )  Deflection  (deg) 

Cf/C- 

Flaps  (plain)  83.6  40 

Stabilizer  30.0 

up  30° 

Elevator  14.1  down20° 

Aileron  5.4  20° 

Rudder  6.0  15° 

Mass  and  Inertia  Characteristics  for  Data  of  this  Report 

.24 

.23 

.18 

.  391  base 
.  453  tip 

Gross  Weight 

2948  pounds 

Center  of  Gravity 

for  the  flight  tests,  27.4%  mac 

for  the  tabulated  derivatives,  25.  0% 

mac 

I 

X 

1284  slug-ft3 

I 

y 

2773  slug-ft3 

I 

z 

3235  slug-ft3 

Flight  Characteristics 

Condition  I 

Condition  II 

Density  Altitude 

5000  ft 

5000  ft 

Velocity 

240  ft/  sec  TAS 

144  ft/  sec  TAS 

Flaps 

0° 

20° 

Lift  Coefficient 

.28 

.75 

Thrust  Coefficient,  T  1 

'  (est. )  .  023 

.  057 
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Lateral -Directional  Degrees  of  Freedom 
1 .  Flight  test  derivatives 

The  matching,  on  an  analog  computer,  of  the  flight -recorded  transient 
responses  has  resulted  in  the  derivatives  tabulated  in  column  7  of  Table  2. 

The  basis  for  the  analog  computer  model  was  the  set  of  more-or-less 
standard  equations  below: 


;  -  Yv)  40 

+  Ar 

v  Ac? 

=  —  Y.  A  6 

V  5r  r 

-L(S40 

-  L  Ar  + 

r 

s(s  -  L  )A;+ 
P 

=  L.  A  6  +  L.  A  6 

6r  r  Oa  a 

-NfJ40 

*  (s  -  N  )Ar  - 
r 

N  s  Aop 

P 

=  N.  A  6  +  N.  A  5 

Or  r  oa  a 

The  various  nondimen sional  derivatives  are  found,  of  course,  from 
the  parameters  of  these  equations  by  applying  the  dimensionalizing  factors: 
airplane  mass,  moment  of  inertia,  wing  area,  span,  dynamic  pressure,  etc.  . 
The  moments  of  inertia  have  been  separately  determined  by  oscillating  the 
airplane  in  different  pendulum  modes,  with  the  following  results 

I  =  1284  slug -ft8 
xx 

1  =  3235  slug -ft2 

zz  6 

These  are  the  corrected  values  for  the  weight  and  loading  conditions 
of  the  flight  tests.  The  product  of  inertia,  I  ,  which  is  omitted  from  the 
above,  is  negligible  for  this  airplane  in  these  flight  conditions  where  the  X 
principal  axis  almost  coincides  with  the  flight  path. 

A  few  brief  comments  about  the  analog  matches  are  in  order  (Figures  2 
and  3).  In  the  first  place,  two  of  the  derivatives  were  determined  from  souices 
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outside  the  flight  test  data.  They  are  C^^  and  Cnga  evaluated  from  the  wind- 
tunnel  data  of  Reference  3,  It  was  not  feasible,  because  of  the  short  roll  time- 
constant  of  the  airplane,  to  distinguish  in  the  analog  matching  between  and 
and  between  Cnp  and  Cng&.  Because  of  this  redundancy  in  the  analog 
matches,  over  fairly  wide  ranges  of  values,  the  effects  of  changes  in  the  aile¬ 
ron  derivatives  would  be  as  follows 


%~C*6a 
C/j3  ~  C^6a 
C^r  ~  C^6a 


C; 

*'P 

C-t6a 


AC 


It  is  considered,  however,  that  the  aileron  control  derivatives  would  be 
the  most  reliable  of  any  determined  from  the  wind-tunnel  data,  and  that  the  pro-  • 
cedure  followed  is  therefore  the  best  that  can  be  done. 

The  roll -due -to -yaw  derivative,  CiY>  has  been  determined  from  a  re¬ 
quirement  to  match  the  spiral  mode  stability  observed  in  special  flight  tests 
for  the  purpose.  In  Condition  I,  the  spiral  mode  was  exactly  neutral;  in  Condi¬ 
tion  II,  the  mode  was  unstable,  Ts  =  10  secs.  These  characteristics  are  matched 
by  the  derivatives  tabulated.  In  the  transients  of  Figures  2  and  3  the  mode  is 
not  excited  and  not  visible. 

In  addition  to  the  transient  responses  and  the  special  spiral  stability 
runs,  steady  sideslips  were  performed  to  provide  equilibrium  equations  in¬ 
volving  static  stability  and  control  derivatives.  Although  in  principle  these 
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data  are  capable  of  removing  the  redundancy  in  the  determination  of  deriva¬ 
tives  from  flight  transients;  in  practice,  the  necessary  solutions  were  poorly 
conditioned,  and  the  sideslip  data  could  only  be  used  to  check  for  consistency 
the  derivatives  determined  from  the  matching  of  transients.  This  check  was 
made,  and  indeed  the  steady  sideslip  equilibrium  equations  are  satisfied  to 
within  the  accuracies  of  derivatives  estimated  below. 

The  Individual  Derivatives  (Table  2) 

The  stability  and  control  derivatives  tabulated  are  discussed  individually 

below. 

1)  The  aileron  effectiveness  has  been  determined  for  Condition  I 

directly  from  full-scale  wind-tunnel  data.  In  the  table,  the  asterisk 
indicates  that  the  flight  value  is  not  an  independent  result. 

For  Condition  II,  wind-tunnel  data  are  not  available  and  the 
flight  value  is  arbitrary.  Both  C^a  and  C.^  should,  however,  be 
little  affected  by  flap  deflection;  and  it  is  considered  to  be  a  con¬ 
firming  factor  that,  starting  with  the  same  Ct^a  in  both  Condi¬ 
tions  I  and  II,  essentially  the  same  values  of  C^  are  found. 

The  control  derivative  is  not  involved  in  the  estimation  of  loads, 
due  to  turbulence,  and  so  its  value  has  not  been  estimated  by  the 
various  other  references. 
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2)  Cn^a  The  aileron  yaw,  for  Condition  I,  has  been  evaluated  from  the 

wind-tunnel  data.  Again,  the  asterisk  indicates  that  the  flight  value 
is  not  independent. 

For  Condition  II,  Cnga  should  be  different,  since  it  would  be 
affected  by  flap  deflection  and  difference  of  lift  coefficient.  As  dis¬ 
cussed  above,  it  can  not,  however,  in  matching  transients,  be  dis¬ 
tinguished  from  Cnp>  and  what  can  be  determined  is  only  a  linear 
combination  of  the  two.  The  particular  values  tabulated  for  flight 
produce  a  proper  match,  but  the  asterisk  indicates  that  they  have 
not  been  separately  identified. 

3)  Cn^r  The  rudder  effectiveness  is  determined  uniquely  in  the  match 

of  flight  responses  by  the  size  of  yaw  rate  response  to  the  rudder 
pulse.  The  first  three  peaks  in  Figures  2  and  3  have  been  clipped 
by  saturation  of  an  instrumentation  amplifier,  and  so  they  are  ig¬ 
nored  in  making  the  match.  The  resulting  accuracy  is  considered 
to  be  about  10  percent. 

This  derivative  is  not  involved  in  the  loads  due  to  turbulence, 
and  so  estimates  by  the  various  references  are  omitted.  The  wind- 
tunnel  value  is  available,  however,  for  Condition  I,  and  the  agree¬ 
ment  is  excellent. 

4)  The  roll  due  to  rudder,  over  a  reasonable  range  of  values,  has 
been  found  to  affect  the  transient  responses  very  little  on  the  Navion. 
Although  a  simple  estimate  was  made  and  provided  on  the  analog 
model,  no  real  determination  was  feasible  and  no  estimate  of  ac¬ 
curacy  is  possible.  For  these  reasons  the  derivative  is  omitted 
entirely  from  the  table. 
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5)  Ci 


6) 


7)  C,, 


8)  Cn 


Although  the  roll  damping  is  not  explicitly  involved  in  the  simpli¬ 
fied  method  of  loads  estimation,  its  value  is  determined  quite  ac¬ 
curately  (say  5  percent,  relative  to  C^^a)  by  the  match  of  "steady- 
state"  roll  rate  response  to  the  aileron  doublet. 

The  values  of  are  adequately  predicted  by  the  various  refer - 
P 

ences,  except  that  Datcom  (Reference  4)  underestimates  by  about 
3  0  percent. 

Over  the  range  of  plausible  values  for  the  Navion,  this  roll -due  - 
to-yaw  derivative  has  no  appreciable  affect  on  the  transient  responses 
of  Figures  2  and  3.  The  values  tabulated  under  flight  have  been  calcu¬ 
lated  so  as  to  match  the  spiral  stability  observed  in  special  flight  tests 
for  the  purpose.  The  accuracy  is  perhaps  20  percent. 

At  least  for  Condition  I,  the  various  estimates  appear  satisfactory, 
except  that  Datcom  (Reference  4)  overestimates  considerably, 
j  The  dihedral  effect  is  evaluated,  under  flight,  by  matching  the 
size  of  the  roll-rate  response  to  the  rudder  pulse.  This  effectively 
controls  also  the  Dutch-Roll  mode -shape,  (p/  r)^.  The  accuracy  is 
considered  to  be  quite  good,  say,  15  percent. 

The  estimates  of  the  various  references  are  reasonable,  with 
Datcom  and  Etkin  (References  4  and  6)  overestimating  about  50  per¬ 
cent.  Even  the  wind-tunnel  value  is  too  large! 

2  The  directional  stability  parameter  is  adjusted  to  match  the 
period  of  the  Dutch-Roll  mode.  It  is  well  and  uniquely  determined, 
with  probable  accuracy  within  10  percent. 

There  is  a  wide  variation  of  values  from  the  different  refer¬ 
ences.  Seckel  (Reference  7)  seriously  underestimates  Cn£  >  and 
even  the  wind-tunnel  data  yield  a  value  too  high  by  about  25  percent. 
Inconsistency  of  the  reference  estimates  might,  in  fact,  be  expected 
from  the  warnings  of  Campbell  and  McKinney  (Reference  2)  that  the 
important  tail  contribution  can  only  be  evaluated  by  reference  to 
empirical  data. 


12, 

Although  the  details  are  by  no  means  clear,  it  seems  likely 
that  the  Seckel  formulae  underestimate  the  stable  vertical  tail  con¬ 
tribution  and  overestimate  the  unstable  fuselage  part.  Etkin's 
method  (Reference  6  )  seems  to  overestimate  seriously  the  tail 
contribution.  A  large  part  of  the  problem  may  be  inability  to 
evaluate  the  sidewash  factor,  dff/ d/3 ,  as  indeed  suggested  by 
Campbell  and  McKinney. 

The  Etkin  text  (Reference  6},  which  appears  to  overestimate 
the  tail  contribution,  leaves  to  judgment  the  evaluation  of  vertical 
taU  area  and  aspect  ratio,  and  so  different  users  of  this  reference 
would  find  different  values  of  the  static  stability. 

9)  Cn^  The  yaw  damping  is  adjusted  to  match  the  damping  of  the  Dutch- 
Roll  mode  in  the  flight  transients.  Its  effect  is  strong  and  direct, 
and  the  accuracy  is  considered  to  be  about  15  percent. 

Discrepancies  among  the  references  are  almost  as  bad  as  for 
Cn^ .  Datcom  (Reference  4)  apparently  overestimates  by  about  50 
percent.  The  statements  of  Campbell  and  McKinney  (Reference  2  ) 
are  again  seemingly  justified,  and  again  the  culprit  may  be  side- 
wash  effects.  Etkin  (Reference  6)  again  seems  to  overestimate 
the  tail  contribution. 

10)  Cn  The  adverse  yaw  derivative  is  determined  by  the  size  of  the 

P 

yaw  rate  response  to  the  aileron  doublet.  For  Condition  I,  where 

a  wind-tunnel  value  of  Cn§a  was  used,  the  Cttp  value  is  valid  and 

should  be  accurate  to  within  about  10  percent.  For  Condition  II, 

as  discussed  earlier,  the  individual  values  of  Cn  and  Cnc  are 

p  uoa 

not  really  known,  as  indicated  by  the  asterisk  in  the  table,  al¬ 
though  a  linear  combination  is  determined  and  represented  in  the 
analog  model. 
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11) 


None  of  the  various  references  are  very  precise  on  a  percentage 
basis,  although  in  the  range  of  values  tabulated,  the  errors  have  little 
influence  on  the  airplane  behavior.  Perkins  and  Hage  (Reference  5) 
predicts  the  right  value  by  neglecting  a  drag  contribution! 

The  flight  value  tabulated  for  the  side -force  derivative  is  derived 
from  the  bank-angle  to  sideslip  ratio  in.  steady  sideslips.  The  reduc¬ 
tion  to  Cy^  involves,  of  course,  the  side -force  due  to  rudder,  which 
is  estimated  from  Cngv  and  the  tail  length  for  rudder  forces.  The 
tabulated  value  is  considered  accurate  within,  say,  10  percent.  For 
Condition  II,  sideslip  data  were  not  available,  and  the  value  was  arbi¬ 
trarily  set  the  same  as  for  Condition  L  The  fact  that  no  valid  deter¬ 
mination  was  made  is  indicated  by  the  asterisk  for  that  case. 

No  certain  explanation  of  the  general  discrepancy  between  the 
various  references  and  the  flight  value  is  available,  although  it  may 
be  assumed  that  the  difference  is  due  to  propeller,  fuselage,  and 
wing  contributions  not  properly  considered.  The  effects  on  the  air¬ 
plane  would  not  be  significant. 


Longitudinal  Degrees  of  Freedom 

Longitudinal  stability  and  control  derivatives  estimated  by  the  various 
references,  and  derived  from  wind-tunnel  data  and  flight  test  transients,  are 
shown  in  Table  3.  The  wind-tunnel  data  would  normally  be  the  standard  against 
which  to  compare  the  reference  estimates,  but  in  this  case,  the  special  flight 
test  data,  with  certain  qualifications,  are  considered  to  be  preferable. 

The  flight  test  derivatives  are  derived  from  trial -and -error  analog - 
matching  of  the  transients  shown  in  Figure  4. 
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The  ba«i*»  for  the  analog  computer  model  was  the  following  equations 


's  +  (D  -  T  f  AV  +  CD  -g]An+gAl  =  0 
v  v  '  a 


L  L 

~  lv  +  is  +  ~Hq  -  she  =  -  ~  h6  tb 
V  V  V  e 

o  o  o  c 


-M  AV  •  [M.  s+M  ]Aflr  +  s  [s  -  MX]  A8  =  M,  A  6 
v  a  a  8  o_  e 


V 

Ati  =  ~  (s  A  8  -  s  Aff) 
g 

The  dimensional  derivatives  above  are  converted,  of  course,  to  their 
nondimensional  counterparts  by  applying  the  factors  of  mass,  inertia,  area, 
etc. 

For  the  kind  of  znotion  involved  in  the  flight  tests,  where  A  V  was 
practically  negligible,  it  is  well  known  *hat  the  parameters,  which  are  in 
principle  unique,  are 


Cl^  -  slope  of  lift -curve 


C™.  -  elevator  effectiveness 

mg 


CLaCmq 

+  — — - —  -  short -period  frequency  parameter 

a 


C  +  C 
m .  m 

<*  q 


-  short -period  damping  parameter 
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For  the  Navion  airplane,  besides  the  redundancy  in  the  individual 
derivatives,  there  is  the  added  difficulty  that  the  short  period  mode  is  of 
high  frequency  and  heavy  damping,  making  it  difficult  to  distinguish  in  prac¬ 
tice  between  the  last  three  parameters  above.  Ratios  of  any  two  to  the  third 
are  easily  found,  but  unique  values  of  all  three  are  quite  difficult  to  resolve 
by  simple  matching.  For  this  reason,  Cm^  was  separately  evaluated  by  a 
special  flight  test  involving  a  difference  in  elevator  position  to  trim  at  two 
different  CG  positions.  This  independent  determination  of  Cm^  resolved 
the  above  redundancy  problem  in  a  satisfactory  way.  As  the  Figure  4  shows, 
the  tabulated  values  provide  an  excellent  match  to  the  flight  test  responses. 

The  analog  model  of  the  longitudinal  dynamics,  of  course,  contains 
the  pitching  moment -of -inertia,  which  has  been  separately  determined  from 
special  oscillation  tests.  The  value,  for  the  flight  condition  and  loading  of 
the  tests,  was 


I  =  2773  slug -ft3 

yy 

Although  only  the  "short -period"  parameters  above  were  varied  in 
the  trial  and  error  matching,  all  three  longitudinal  degrees  of  freedom  were 
provided  in  the  model.  The  parameters  of  the  drag  equation  and  the  coeffi¬ 
cients  of  velocity  change,  A  V  ,  were  pure  estimates  and  no  attempt  was 
made  to  check  their  values.  This  detail  is  not  significant,  since  it  could 
be  seen  in  the  data  that  the  velocity  changes,  during  the  transients  shown, 
were  very  small  and  seemingly  negligible  in  any  case. 

Discussion  of  the  individual  derivatives  follows. 

1)  Cm«.  The  elevator  effectiveness  controls,  in  principle,  the  size  of 

load  -factor  and  pitch -rate  responses.  However,  the  static  stability 


parameter 
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'm 


at 


C.  C 
L  m 

a  q 

4}j 


is  also  involved,  and  in  the  flight  test  transients,  their  effects  can 

not  really  be  separated,  as  previously  noted.  This  was  resolved  by 

a  separate  flight  test  for  Cm,  ,  with  the  result  tabulated.  The  ac - 

*e 

curacy  of  this  special  result  is  perhaps  10  percent,  and  so  the  per¬ 
fect  agreement  with  the  wind  tunnel  value  is  fortuitous. 

The  control  derivative  is  not  required  in  the  estimation  of  loads 
due  to  turbulence,  and  so  values  from  the  various  references  are  not 
tabulated. 


2)  Ci  The  important  slope  of  the  lift -curve  is  determined  by  matching 
-'a 

the  size  of  the  load -factor  responses  in  Figure  4.  The  accuracy, 
relative  to  the  control  effectiveness,  should  be  excellent,  perhaps 
10  percent  or  better. 

The  estimates  from  the  reference  textbooks  may  be  low  because 
they  only  count  wing  lift.  But  of  course  for  wing  loads,  one  should 
only  count  the  wing  contribution,  r  nd  so  for  that  purpose  the  esti¬ 
mates  may  be  satisfactory. 

C  t  C 
L  m 

(X  Q 

3)  Cm/V  +  - ^ — —  The  static  stability,  or  frequency,  parameter  could 


'mQf 


ideally  be  determined,  in  the  match,  by  the  timing  of  successive  peaks 
in  the  pitch  rate  responses.  Especially  in  Condition  I,  the  natural  fre¬ 
quency  is  so  high  that  the  method  is  insensitive;  and  so  the  parameter 
is  used  to  control  the  size  of  pitch  rate  response  instead,  for  the  given 


m6  ’ 
°e 

cent. 


Done  in  this  way,  the  accuracy  should  be  excellent  -  say  10  per- 
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For  Condition  II,  where  the  natural  frequency  was  lower,  a  partly 
independent  determination  could  be  made,  although  not  with  good  ac¬ 
curacy.  In  this  case  Cm^  was  set  to  the  wind-tunnel  value,  which  was 
validated  in  Condition  I,  and  the  resulting  accuracy  of  the  static  sta¬ 
bility  parameter  is  estimated  at  5  or  10  percent.  A  confirming  feature 
of  the  matched  responses  is  the  perfect  timing  of  the  response  peaks  in 
both  acceleration  and  pitch  rate. 


The  various  references  exhibit  an  unpleasantly  large  variation  in 
their  estimates  of  this  important  stability  parameter.  The  largest 
error,  by  Seckel  (Reference  7),  amounts  to  an  equivalent  CG  shift  of 
12  percent  mac;  while  the  smallest  error,  by  Datcom  (Reference  4) 
represents  about  2  percent  mac.  It  seems  probable  that  the  large 
errors  are  due  to  problems  of  stimating  the  important  stabilizer 
effectiveness,  Cmjt.  This  would  be  significant  in  the  estimation  of 
horizontal  tail  loads  due  to  turbulence,  as  for  example  in  Reference  1. 

Under  wind-tunnel,  the  overall  value  is  not  directly  available.  It 
has  been  estimated  using  only  the  tail  contribution  to  CmQ,  from  the 
formula 


A,  C 


mr 


=  2  —  C 


mit 


4)  C  +  C 
m .  m 
a  q 


This  qualification  is  indicated  by  asterisks. 

The  rate  damping  parameter  is  determined  in  matching 

by  ratios  between  peak  amplitudes  in  the  pitch  rate  responses, 
Figure  4.  They  are  quite  sensitive  to  the  parameter,  ahd  so  the 
accuracy  should  be  good,  say,  1  0  percent. 

Although  the  textbook  formulas  predict  too  small  values,  the 
Datcom  estimates  are  a  bit  too  large;  but  in  fact  all  the  estimates 
are  probably  satisfactory. 
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CONCLUSIONS 


Based  on  the  study  of  stability  derivatives  of  the  Navion,  as  estimated 
by  various  references  and  compared  to  full-scale  wind-tunnel  and  flight  test 
results,  the  following  conclusions  are  drawn.  They  are  presumed  to  apply 
generally  to  the  estimation  of  those  parameters  for  light  single -engine  pro¬ 
peller  airplanes  of  conventional  configuration  with  high  aspect  ratio  straight 
wing  and  tail,  to  flight  conditions  at  low  angle  of  attack  and  where  power  and 
slipstream  effects  are  small. 

1)  The  important  airplane  slope  of  the  lift -curve,  Cl^,  is  under¬ 
estimated  by  the  textbook  references,  probably  because  of  neglecting  body 


and  tail  lift.  They  may,  however,  be  satisfactory  for  estimating  wing  loads. 

2)  The  static  stability  in  angle -of -attack,  Cm^  + - ^ — —  varies 

considerably  as  estimated  by  the  different  references.  In  terms  of  equivalent 


CG  position,  the  errors  vary  from  2  percent  to  12  percent  chord.  The  dif¬ 
ficulty  appears  to  be  with  the  important  stabilizer  effectiveness. 

3)  Tne  rate  damping  term,  Cm  .  +  Cm  ,  appears  to  be  underestimated 

a  q 

by  all  the  textbook  references  -  by  about  20  percent. 

4)  The  directional  stability,  Cnp  ,  is  predicted  well  by  the  references 
except  one,  where  an  overestimate  of  the  fuselage  contribution  and  a  somewhat 
low  tail  term  combine  to  produce  a  very  low  overall  value. 

5)  The  yaw  damping,  Cnr,  is  overestimated  badly  (~  50  percent)  by  two 
of  the  references.  The  difficulties  with  Cnr  and  Cn^  are  probably  at  least  partly 
due  to  sidewash  effects  on  the  flow  at  the  vertical  tail.  These  are  not  adequately 
treated  anywhere,  and  the  statements  of  Campbell  and  McKinney,  about  the  need 
for  test  data  on  a  similar  configuration,  seem  to  be  confirmed. 

6)  The  roll  damping,  ,  is  adequately  predicted  by  the  references 

P 

save  one,  which  comes  in  low  by  about  3  0  percent. 

7)  The  dihedral  effect,  ,  is  overestimated  about  50  percent  by  two 
of  the  references.  The  others  are  satisfactory. 
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8)  The  roll  due  to  yaw,  C/^,  is  adequately  predicted  by  the  references 
except  one  which  overestimates  it  by  some  90  percent. 

9)  The  yaw  due  to  roll,  Cnp,  exhibits  large  fractional  errors  by  the 
references;  but  in  the  flight  condition  considered,  the  level  of  the  derivative 
is  small,  so  that  the  errors  are  not  very  significant. 

10)  The  side  force  derivative,  Cy^ ,  is  grossly  underestimated  by  all 
the  references;  probably,  however,  without  much  affect  on  the  airplane, 

11)  Particular  values  of  all  the  derivatives  cited  above  are  listed  in 
Tables  2  and  3.  Estimates  judged  to  be  unacceptably  erroneous  are  circled 
for  easy  identification. 

12)  It  is  amply  clear  that  the  accuracy  of  any  of  the  cited  references, 
followed  literally,  would  be  inferior  to  that  of  an  "old  hand"  at  airplane  aero¬ 
dynamics,  who  would  be  guided  by  experience  and  related  test  data  to  modify 
the  predictions  in  "ad-hoc"  fashion. 
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TABLE  II  -  LATERAL  DIRECTIONAL  DYNAMIC  DERIVATIVES 


Derivative 

Textbooks 

NACA 

TR  1098 

USAF 

Datcom 

Wind 

Tunnel 

Flight 

Test 

Perkins 
and  Hage 

Etkin 

Seckel 

cy/3 

- 

-.28 

-.20 

-.20 

- 

-.35 

77 

-.  61 

-.61* 

(Cy  _  ) 

yP  Vt 

- 

- 

-.  18 

-.18 

- 

- 

- 

- 

062 

-.062 

-.  077 

-.077 

-.  072 

-.  100  | 

-.109 

-.097 

-.  067 

-.  051 

| 

C^r 

.070 

.  195 

.084 

.  067 

.  141 

.062 

■ggym 

- 

.069 

.27 

% 

-.46 

-.  46 

-.39 

-.42 

-.42 

-.42 

-.30  j 

- 

-.46 

-.48 

S 

.  073 

.073 

.094 

BhH 

- 

.  082 

.  109 

.  109 

.  086 

.  084 

(<Vvt 

.  109 

.  109 

.139 

- 

- 

- 

cn 

nr 

-.087 

-.096 

>141 J 

-.065 

-.  076 

- 

120 1- ■ 

- 

-.  088 

-.  163 

'Wvt 

084 

-.  084 

-.  139 

-.  061 

-.  061 

- 

- 

- 

- 

S 

-.  035 

-.098 

-.017 

-.  010 

-.030 

- 

C^6a 

- 

- 

- 

- 

.  152 

.  152=:= 

.iso* 

Cn6a 

- 

- 

- 

- 

-.0047 

Cn6r 

- 

- 

- 

- 

- 

-.  077 

-.  075 

-.  093 

NOTES 

1,  In  each  box,  the  number  in  the  upper  left  corner  is  for  Condition  I  and  the 
number  in  the  lower  right  corner  is  for  Condition  II. 

2,  Theenboxed  numbers  are  estimates  which  are  considered  to  be  unacceptably 
erroneous, 

3,  The  asterisks  signify  numbers  which  were  not  found  independently  during  the 
analog  matching  procedure  (see  Discussion). 
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TABLE  III  -  LONGITUDINAL  DYNAMIC  DERIVATIVES 


Derivative 

Textbooks 

USAF 

Datcom 

Wind 

Tunnel 

1 

Perkins 
and  Hage 

Etkin 

Seckel 

Flight 

Test 

4.36 

4.25 

4,54 

5.50 

4.52 

6.  04 

6.40 

°m  a 

-.83 

-.715 

-.872 

-.545 

-.545 

-1.24 

-.95 

-1.  03 

- 

Cmd 

-3.0 

-4.91 

-5.23 

-4.91 

-4.91 

-6.58 

- 

Cmq 

-9.6 

-9.  75 

-9.  75 

-9.50 

-9.50 

-13.29 

- 

- 

Cm5e 

- 

- 

- 

- 

-1.42 

-1. 55 

-1.42 

-1.55* 

Cmit 

-1. 82 

-1.68 

-1.68 

-1.  72 

-1.72 

wm 

1 

- 

CL*CMq 

4U 

-1.08  J 

ESalB 

mamA‘ 

-1.68 

-1.27* 

-  1.35* 

-1.55 

-1.37 

■0 

cm<$  +  Cmq 

-12.6 

-14.7 

-14.4 

-14.4 

-19.  8 

- 

-18.3 

-15.5 

NOTES 


1.  In  each  box,  the  number  in  the  upper  left  corner  is  for  Condition  I 
and  the  number  in  the  lower  right  corner  is  for  Condition  II. 

2.  Theenboxed  numbers  are  estimates  which  are  considered  to  be 
unacceptably  erroneous. 

3.  The  asterisks  signify  numbers  which  were  not  found  independently 
(see  Discussion). 
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riGURE  2  ANALOG  COMPUTER  MATCH  OF  CONTROL  RESPONSES;  CONDITION  I 

FLAPS  UP  ,  CL«.28  ’ 

LATERAL  -  DIRECTIONAL  MOTION 
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FIGURE  3  ANALOG  COMPUTER  MATCH  OF  CONTROL  RESPONSES  ;  CONDITION  IE , 

8f  *  20°  ,  CL  »  .75 
LATERAL  -  DIRECTIONAL  MOTION 
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FIGURE  4  ANALOG  COMPUTER  MATCH  OF  CONTROL  RESPONSE 

LONGITUDINAL  MOTION 


